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Motivated by the recent discovery of two Q ~ collective modes [Y. Li et al, Nature 468, 283 (2010); 
ibid Nat. Phys. 8, 404 (2012)] in the single-layer HgBa2Cu04+<5, which are often taken as evidence of the 
orbital current origin of a pseudogap, we examine an alternative and assumption-free scenario constrained by 
first-principle calculations. We find that in addition to the common Cu02 band, a hybridized Hg-0 state is 
present in the vicinity of the Fermi level, and that it contributes to the low-energy ground state of this system. 
We calculate the spin-excitation spectrum based on the random-phase- approximation in the superconducting 
state using a two band model and show that a collective mode in the multi-orbital channel arises at Q = 0. 
This mode splits in energy, yet remains at Q ~ as the pseudogap develops breaking both translational and 
time-reversal symmetries in the spin-channel, and thus gap out the electrnic states. The observations of the 
dynamical mode and static moment in the pseudogap state are in good accord with experimental observations. 
A verifiable prediction of this proposal is the development of magnetic order also in the Hg-0 layer. Detection 
of Hg-0 band via optical study, or magnetic moment in the Hg-0 layer will be tests of this calculation. 

PACS numbers: 74.72.Gh,74.72.Kf,75.25.-j 



Unraveling the nature and mechanism of the pseudogap has 
remained a steady theme of the research of cuprates, 1 and 
has recently been extended to pnictide 2 and heavy fermion 
systems. 3 Within cuprates, the evidence of pseudogaps has 
emerged to be highly contradictory in different experimen- 
tal probes as well as in different materials. While the orig- 
inal idea of the 'pseudogap' stems from the observations of 
anomalous (non mean-field like) temperature dependencies in 
various bulk measurements, 4 and incoherent spectral weight 
properties 5 in spectroscopies that set in below a characteristic 
temperature T* 9 recent magneto-transport measurements have 
indicated the presence of some form of density wave origin of 
it.-. Focusing on the broken symmetry state of the pseudogap 
from various inelastic neutron scattering (INS) data, its nature 
apparently seems to be strongly material dependent within the 
cuprate family as follows. (1) Electron doped cuprates consis- 
tently demonstrate the presence of a commensurate antiferro- 
magnetic order up to the superconducting (SC) region. 7 (2) In 
single layer hole doped La2_ ;E (Ba/Sr) ;E Cu04 (LB /SCO), the 
commensurate order is observed to become incommensurate 
with doping (with a so-called 'hour-glass' dispersion in the 
spin-excitation spectrum), and rotates along the Cu-0 bond 
direction as the SC dome develops. 8 Such phenomena have 
been taken as evidence of the 'stripe' order origin of the pseu- 
dogap in these systems.^ (3) In YBa 2 Cu 3 06+5 (YBCO), the 
incommensurate spin-excitation further exhibits an in-plane 
anisotropy, having stronger intensity along the a-bond direc- 
tion than the 6-direction, a fact which has been interpreted as 
the emergence of electronic nematic order. 10 (4) In addition 
to the 'hour-glass' pattern, which is a trademark of all these 
hole doped systems, a Q ~ mode is observed in single layer 
HgBa2Cu04 + £ (Hg 1201), 11 which further splits in energy.— 
Such an Ising-like mode can be expected within a circulating 
orbital current model.— 

Understanding such diverse material dependence of the 
spin-excitation properties in cuprates within a single model 
has remained a challenge. However, some efforts to obtain a 



unified model interpretation should be mentioned. The exten- 
sion of the model of coexistence of an (either long- or short- 
ranged) spin-density wave (SDW) and d-wave superconduc- 
tivity for electron doped cuprates to the hole doped side is one 
of these successful approaches. Within this model, it has 
been shown that the commensurate mode is lifted to higher 
energy by the SC gap with a downward dispersion reaching to 
zero energy due to the nodal d-wave momentum (fe) depen- 
dence of the SC gap.i 5 ^ In other words, the spin-excitation 
manifests itself as static incommensurate peaks via coupling 
to the nodal d-wave SC gap. This model reproduces the 'hour- 
glass' pattern, a 45° rotation of the incommensurate spec- 
tra in all these systems. Furthermore, the same model has 
been applied to quantitatively explain the additional in-plane 
anisotropy in YBCO by incorporating the inter-layer coupling 
between the Cu02 plane and the metallic unidirectional CuO 
chain plane which is present at all finite dopings. 17 This cal- 
culation does not require any 'spontaneous' nematic order. In 
the present paper, we extend the model to the Hgl201 system 
to explain the coexistence of both the 'hour-glass' pattern and 
the Q ~ modes. 

Looking into the first-principles band structure,— - — we find 
that there exists a Hg-0 band which lies so close to the Fermi 
level (Ep) that it can contribute to the low-energy physics in 
this compound. Constrained by this band- structure property, 
we present an assumption-free two band tight-binding model 
including the Cu02 plane and Hg-atomic plane, see Fig. [Ha). 
We compute the spin excitation spectrum to show that due to 
the interaction between these two plane states, a Q ~ col- 
lective mode develops in the multi-orbital spin channel. We 
further demonstrate that when a SDW order or any other sim- 
ilar translational symmetry breaking order sets in the CuO 2 
plane in the pseudogap regime, it induces a magnetic moment 
to the Hg-0 state which splits the Q ~ mode in energy, but 
not in momentum. The observations of the dynamical mode 
and associated static moment in the pseudogap state are in 
good accord with experimental observations . n i 12 i 21 
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FIG. 1 . (Color online) (a) Crystal structure and tight-binding hopping 
parameters of the Hgl201 system, (b) Non-interacting tight-binding 
FS (red solid line) plotted on top of the first-principles result.—— (c) 
Corresponding dispersion plotted along high- symmetry directions in 
solid colors and compared with ab-initio dispersions. 
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FIG. 2. (Color online) Imaginary part of the BCS-RPA susceptibil- 
ity in the paramagnetic state plotted along the diagonal direction as a 
function of excitation energy, (a) The intra-atomic component for the 
Cu02 band reveals strong peak in the Q2 ~ (7r, tt) region. Although 
in the bare level, this component does not contain any other intensity 
except at Q2, its RPA value involves weak features near Qi ~ 
due to mixing with other terms within the tensor form of the RPA 
denominator (discussed in the text), (b) The inter- atomic RPA sus- 
ceptibility between the CUO2 and Hg-0 bands shows intensity at Q±. 
The zone periodicity of the intensity between Q± and Q2 even in the 
paramagnetic state strongly suggests that the Qi mode is unstable to 
a (fluctuating) magnetic ground state with a modulation of Q2. (c) 
The total RPA susceptibility yields strong intensity at both Q vectors 
of different origins, in agreement with experimental data.— 



Tight-binding Hamiltonian:- Unlike in other cuprates 
where the Cu-O-0 layer(s) dominate the most interesting low- 
energy electronic states, in the Hg-based cuprates the hy- 
bridization of the reservoir atom Hg with the apical oxygens 
and Ba is strong, making the Hg-0 derived band cross or 
stay close to Ep. Both x-ray photoemission data^ and first- 
principles^ - — calculations have demonstrated the presence of 
the bottom of the Hg-0 band lying as close as 0.09-0.13 eV 
above Ep at the 'X' point which moves and even crosses 
Ep with doping, interaction, or with an increasing number 
of Q1O2 layers. The first-principles calculations also indi- 
cate that individual Q1O2 and Hg-0 bands are highly two- 
dimensional, and the inter-layer coupling between them is 
prominent. Based on these considerations, we derive a re- 
alistic two-band tight-binding model for Hgl201, in which 
the parameters are deduced by fitting to first-principle disper- 
sion, without any adjustments. The obtained non-interacting 
Hamiltonian is 

Hh= (t*t™\. (1) 

The dispersions £ifc and ^2k are for the Cu02 and Hg-0 
states, respectively, where £i2fc is the inter-layer coupling be- 
tween them, see Fig. [Ha), in which ^ = —2ti((j) x + <j) y ) — 

^fafa - 2% (fax + fay) ~ ^'"(faxfa + fay fa) ~ f*i, and 

fi2fc = - 2ti2 (f> z /2 , where fa x / y / z = cos(ak x/y/z ). The 
corresponding values of the tight-binding parameters U and 
chemical potentials fii are given in Ref . |23l 

The tight-binding fittings to the first-principles dispersion 
and Fermi surface (FS) are given in Figs. [Tib) andQIc). The 



CUO2 antibonding state is clearly visible , and it only consti- 
tutes the FS. An essential difference between the Cu02 and 
Hg-0 layers is that despite the absence of O-atoms in the Hg- 
plane, see Fig. Uta), first-principles calculations exhibit that 
the latter state (the green line above Ep) originates from the 
hybridization between the Hg and O atoms. This crucial in- 
formation emphasizes that the inter-layer tunneling matrix- 
element t\2 is important and cannot be neglected for any real- 
istic computation. Based on this observation, it is justifiable to 
assume that the inter-layer interaction, V, is also strong, and 
it is the prime term for the development of the Q ~ mode 
in our model. 

Before including magnetic order, we first study the phys- 
ical origin of the Qi ~ and Q2 ~ (tTjTt) modes in the 
paramagnetic state. The widely believed origin of the spin- 
resonance mode at Q 2 is due to the the sign-reversal of the 
d— wave pairing at the magnetic 'hot-spot' on the FS of the 
Q1O2 state Here we show that the Qi mode is also a 
collective phenomenon of different origin that develops in the 
particle-hole continuum of the the inter-orbital channel within 
the random-phase approximation (RPA). Due to the lack of a 
FS in the Hg-0 band, we assume this band to be non-SC in 
this single layer Hgl201 case. Including all these realistic ef- 
fects, we evaluate the orbital dependent non-interacting BCS 
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susceptibility as 

(xo) l jj'(q,Pm) -- 



k,n,vv' 



Gl{Lo n )G v k+q {Lo n +p m ) - F^ n )F v _ k _ q {Lo n + Pm ) 

(2) 

Here G v and F v are the normal and anomalous Green's func- 
tions for the quasiparticle band E% = \/( e fc) 2 + (^50 2 where 
e v h are the eigenvalues of the Hamiltonian given in Eq.Q] If the 
corresponding eigenstates are denoted by ip" (fe) for the i th 
orbital, then the orbital overlap matrix-element can be writ- 



ten as Mn'j'j 



(k, q) = ^\k)^{k + q)^',\k + qW 3 (k). 



ou n and p m are the fermionic and bosonic Matsubara fre- 
quencies, respectively. The SC gap functions are taken as 
A* = A (</) x - (j> y ) and A| = for the Cu0 2 and Hg-0 
bands, respectively, with Ao=31 meV. Finally, we employ a 
RPA correction to obtain the many-body orbital spin suscep- 
tibility as x(Q,v) = Xo(<7,^)[l - UsXoiq^)]' 1 , where U s 
is the interaction tensor in the transverse spin-flip channel de- 
fined in the orbital basis (see Ref. l24j) . 

The individual and the total components of the RPA-BCS 
susceptibility are presented in the SC state in Fig. [2] The 
imaginary part of the intra-orbital component xli f° r the 
CuC>2 state clearly exhibits the lower-branch of the so-called 
'hour-glass' pattern at Q 2 , which is the trademark feature 
of the d-wave SC gap as ubiquitously measured by INS in 
all hole-doped cuprates , 9 i 1Q i 25 including the present Hgl201 
system. 12 In the inter-orbital channel x\h a strong resonance 
mode appears at Qi ~ 0, with a somewhat upward dispers- 
ing branch with vanishing intensity. Tracking down to the 
band- structure details in Fig. Q]c), we reveal that this mode 
originates from the direct excitation gap in the van-Hove sin- 
gularity region (the RPA correction shifts the mode energy to a 
lower value). It is interesting to note that even in the absence 
of any magnetic order, the Qi mode reveals zone periodic- 
ity at (tt, 7r), suggesting that the Hg-atoms tend to magneti- 
cally order by the same wave vector. Since a Qi ~ mode 
develops in the paramagnetic ground state, it is expected to 
survive in the overdoped region. In this context, we recall a 
recent observation of an unusual Raman mode 26 in overdoped 
Hgl201 which can be taken as the persistence of Qi ~ 
above the pseudogap region. Optical absorption spectroscopy 
which measures the direct gap can also be used to test our 
proposal. 

Magnetic ground state:- Next we focus on how a mag- 
netic order can split the Qi ~ mode. Earlier nuclear 
magnetic resonance (NMR) data on the same sample have 
demonstrated the opening of the pseudogap at T* as in the 
INS measurements. 27 Recent neutron diffraction data on the 
same sample 2 ^ as well as in YBCO 28 29 establish that there 
exists a static magnetic order which vanishes above T*, in ad- 
dition to the dynamical mode that saturates above T*. Re- 
lating the static magnetic order to that that which renders 
the FS reconstruction, we work in an (in-plane) double unit 
cell magnetic Brillouin zone connected by the commensu- 
rate SDW wave vector Q 2 = (tt, 7r). If no other symme- 
try is broken, the new magnetic zone is the same for both 
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FIG. 3. (Color online) (a) Single particle spectral weight map in the 
SDW state. The gap opening occurs in C11O2 at Ef at the antinodal 
region, while the smaller magnetic gap in the Hg-0 state commences 
above Ef. Inset: As superconductivity is turned on in the CUO2 
state, the corresponding gap size increases, and for these particular 
gap values, the upper CUO2 magnetic band and lower Hg-0 magnetic 
band becomes close to each other, (b) The details of the energy level 
and the corresponding spin excitation transitions are illustrated in the 
paramagnetic and magnetic states. At Q2, the magnetic transition 
across Ef also involves a momentum transfer, as expected, which is 
not explicitly illustrated in this schematic diagram. 



the Cu02 plane as well as the Hg-plane in the tetragonal 
lattice, and thus, a magnetic moment is induced to the lat- 
ter state. Using the standard Nambu notation, we define 
the two atomic eigenstates in the magnetic zone as ^ = 






lfcf °2fcf C i(fc+Q 2 H' °2(fc+Q 2 )|J' where c ika (c ika ) cre- 
ates (annihilates) an electron with momentum k and spin a 
on the i th atom. In this notation the Hamiltonian presented in 
Eq.[[]becomes a 4x4 matrix H = J2 k ^k H k^k- 



( 



\ 



ilk i\2k -U\m\ \ 

i\2k 6fc -U 2 m 2 

-U 1 m 1 £l(fc+Q 2 ) £l2(fc+Q 2 ) 

-U 2 m 2 £l2(fc+Q 2 ) 6(fc+Q 2 ) / 



•(3) 



Here the order parameters rrii represent the staggered mag- 
netic moments, which in the mean-field level are evalu- 
ated self-consistently as m, = £) fco . a ( C l(fc+Q 2 )<r c ^) = 

T,j^i,k* JT ^t aA ij a > ^)/M- Here the spin-resolved 
spectral function (fc, cr, uj) = —ImGij (fc, cr, uj)/7t, where 
/ is the Fermi function, and cr = a = ±. Self-consistent 
values of the magnetic moments are 0.05 /ig and 0.04/i# in 
the CuC>2 and Hg-0 bands, respectively, and their total value 
is close to the experimental result of 0.1 hb for the Hg 1201 21 
and YBCO 28 ' 29 systems. 

The single-particle spectral function ^2 io . Aa(k, cr, u), 
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FIG. 4. (Color online) (a), (b) Theoretical spin-excitation spectra in 
the coexistence of the SDW and SC ground state, plotted along the 
diagonal, and the (100) directions respectively. The symbols give the 
corresponding experimental data for Hgl201 near optimal doping. 11 
The splitting of the Qi mode in the pseudogap state is more evident 
in (b). (c), (d) Intensity plots of the same experimental results near 
Qi and Q2, respectively. It is evident from the experimental data 
that despite the presence of a tail of the weakly dispersing branch 
in (c), the strong intensity is concentrated near Q\ as in our theo- 
retical spectra in (b). (e) The same theoretical result plotted at Q± 
and Q2, demonstrating the relative intensities of the corresponding 
peaks. Here the two split modes at Q± with equivalent intensities are 
clearly visible. 



plotted along the high- symmetry directions in Fig. Ha) ex- 
poses the nature of gap opening in both bands. The direct 
gap values between t he spin- split states i s determined by U mi 
which increases to y (Umi) 2 + (AjJ 2 in the SC state in the 
Q1O2 layer. (There is no SC gap in the Hg-0 state). In the 
magnetic state, G and F also become 4x4 tensors, and with 
this modification, the calculations of the spin resonance spec- 
trum remain the same as before. Several low-energy magnetic 
transition channels in the particle-hole continuum become ac- 
tive in this case as illustrated in Fig. Ob). Earlier calculations 



have shown that the SDW order in the Q1O2 state gives rise to 
an upward dispersion centering on Q2 which meets the down- 
ward dispersion of the SC origin at Q2 to create the so-called 
'hour-glass' phenomenon. 15 Such a pattern is reproduced here 
in the Q1O2 state as shown in Fig. H^a). Due to the lack of 
a FS in the Hg-0 state, such upward dispersion is absent in 
this band. However, in multi-layered HgBa2Ca n _iCu n 02 n +2 
(with n > 1) the Hg-0 band crosses below Ep, and based on 
the same theoretical argument 15 an additional upward branch 
can be expected there. 

In Fig. 01 we give details of the evolution of two Qi ~ 
modes in the magnetic state. As illustrated in Fig. Ob), the ex- 
citation from the lower-magnetic band of the Q1O2 state to the 
two split Hg-0 bands above Ep creates two spin resonance 
modes. This energy splitting is evident in the spin-excitation 
spectrum plotted along the diagonal and (lOO)-directions in 
Figs.lHa), and|Hb), respectively. We include the experimen- 
tal data (symbols) for the optimally doped Hgl201 sample for 
comparison. 12 The Q\ ~ mode splits by about 20 meV and 
lies at 45 and 55 meV, in accord with experimental values. 
Some discrepancies are clearly visible. The intensity plot of 
the experimental data [reproduced from Ref . Il2l in Figs. BJc), 
and ffid)] reveals a weak dispersion mode with its intensity 
sharply vanishing away from the Q = momentum. We find 
that the tail of the intensity in the theoretical curve disperses 
strongly to higher energy than its experimental counterpart. 
On the other hand, experimental data also demonstrated that 
while the intensity of the Q2 mode vanishes at T* 9 the same 
at Qi presumably merges to the background. This supports 
our postulate that the latter is not directly related to the pseu- 
dogap physics, but arise in the particle-hole continuum of the 
paramagnetic ground state. 



I. CONCLUSION 

In summary, we provide a different and realistic explana- 
tion to the experimental observations of the Ising-like spin ex- 
citation modes in Hgl201 based on the presence of additional 
Hg-0 states close to Ep as established by first-principles and 
photoemission electronic structure considerations. The higher 
energy Qi ~ mode is fairly doping independent, and based 
on our calculation we predict that it should survive to over- 
doping. On the other hand the second mode appears as a re- 
sult of the proximity induced magnetic order in the Hg-atoms 
and hence it is doping and temperature dependent. These 
facts are in accord with experimental data. The similar ob- 
servation of the Qi ~ mode in the YBCO sampl o 28 i 29 can 
also be explained within the same framework due to the pres- 
ence of a metallic chain state in this sample. 17 Furthermore, in 
double layered E^S^CaQ^Og+s (Bi2212), first-principles 
calculations 25 find that the Bi-0 band lies close to the Fermi 
level or may even cross it near the same (tt, 0) point and thus 
can give rise to a similar Q ~ magnetic mode, but its en- 
ergy scale is yet to be explored. However, despite extensive 
neutron studies in other hole-doped cuprates, any signature of 
this Q ~ mode has not yet been reported supporting our 
theory. The gap opening in the Hg-0 state can be tested by 
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scanning tunneling microscopy, optical and x-ray absorption 
spectroscopies, whereas the prediction of magnetic order in 
Hg-0 layer can be verified via NMR, neutron measurements. 
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